Radiative decays: a new flavour filter 
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Abstract 

Radiative decays of the l 3 L>i orbital excitations of the p, to and <fi to the 
scalars /o(1370), /o(1500) and /o(1710) are shown to provide a flavour filter, 
clarifying the extent of glueball mixing in the scalar states. A complementary 
approach to the latter is provided by the radiative decays of the scalar mesons 
to the ground-state vectors p, lo and (f>. Discrimination among different mixing 
scenarios is strong. 
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Introduction 



The observation and confirmation of gluonic degrees of freedom in mesonic states is 
of great significance for QCD, which predicts the existence of glueballs (bound states 
of gluons) and of hybrids (quark- ant iquark-gluon states). Calculations in lattice 
QCD give estimates of the likely masses of glueballs |] and light-quark hybrids @, 
and evidence for the excitation of gluonic degrees of freedom has emerged in several 
processes. Lattice calculations predict that the lightest glueball has J PC = ++ 
and is in the mass range 1.3 to 1.7 GeV. Experimentally J3] there is one more 
++ state in this mass range than can be accomodated by excited quark- ant iquark 
states. The natural inference is that there is a glueball state present [Q, |5|. A clear 
exotic resonance, 7Ti(1600), with quantum numbers J PC = 1 h which cannot be 
accessed by a pure quark- ant iquark state, has been seen || in the r/(958)7r channel 
in the reaction tt~N — > (r/(958)7r)iV and in the p°n~ channel |7|, Q in the reaction 
n~N — > (tt + 7t~7t~)N. Anomalously large hadronic decay modes of light-quark vector 
mesons, that is decays which are predicted to be very small in standard models 
0, are observed in e + e~ annihilation around 1.6 GeV. A favoured explanation 
]TD| , |TT| , m ] is to include vector hybrids which, on the basis of the ^(lGOO) mass, 
are expected to be in this mass region. 

Apart from states with exotic quantum numbers, disentangling hybrids and glueballs 
from quark- ant iquark states is difficult using only hadronic decay channels because 
of mixing. However, radiative transitions offer special opportunities. The coupling 
to the charges and spins of constituents reveals detailed information about wave 
functions and discriminates among models. In the case of gluonic excitations of the 
7r and p, that is hybrids, the spin structure differs from conventional excitations with 
the same overall J . For example, in a hybrid 1 the qq are in a spin-singlet, 
while for the h they are in a spin-triplet: in each case this is the reverse of what 
one is accustomed to. In the case of the scalar mesons a direct measure of their 
electromagnetic couplings gives information about the flavour content of the scalar 
states and could resolve the issue of G-qq mixing. 

We have calculated [13| the rates for the radiative decays V* — *> 7M, where the V* 
are the 25" and ID excitations of the p, u and <fi in the 1.4 to 1.8 GeV mass range 
and the M are positive C-parity qq states in the 1.2 to 1.7 GeV mass range. We 
label the vector states ps, Pd, <^s, ^d, (fis an d 4>d- Some of these vector decays 
are predicted to have branching ratios in excess of 10~ 2 and will be measurable at 
existing and planned facilities, such as e + e~ annihilation by Initial State Radiation 
(ISR) at Babar and Belle, by direct e + e~ annihilation at the upgraded VEPP collider 
at Novosibirsk and by diffractive photoproduction following the upgrade at Jefferson 
Laboratory. 

These radiative transitions serve two purposes, exploring the nature both of the 
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initial excited vector state and that of the resultant meson. In |13| we concentrated 
on the former aspect. Here we emphasise the latter and extend the calculation to 
the radiative decays of the scalars to the ground-state vector mesons. There is again 
great potential for precise measurement of these decays, for example at CLEO via 
the decay J/ip — > 77^ or at COMPASS via central production in high-energy proton 
proton collisions. 

Radiative transitions of p* or uj* to 7M couple directly to the uu±dd content of M. 
In similar vein, the analogous transitions involving 0* couple to the B content of M. 
By comparing the relative rates for isovector (|M) = ^(\uu) - \dd)) and isoscalar 

(|M) = cos# |fi) + sin# -^(\uu) + \dd))) for a set of mesons M of positive C-parity 
and the same J p it is possible in principle to determine the relative amount of 8 and 
non-strange flavours in the M wavefunction and hence to weigh the flavour content 
of the nonet. In the case where a glueball has mixed into the multiplet extending 
it to a decuplet, as is hypothesised @, ^| to be the case for the / (1370), /o(1500) 
and / (1710), such transitions enable the role of the glueball to be disentangled and 
might even be able to estabish the mass of the "quenched" glueball 0. Furthermore, 
as we shall show, radiative decays of the scalar mesons to the ground-state vectors 
p(770) and 0(1020) are also sensitive probes of the glueball mixing in the scalars. 



The model 

The details of the radiative decay calculation are described in [JT3| . Wavefunctions 
for the excited vector mesons are found variationally from the Hamiltonian 



p 2 A a 



H = ^— + ar--— + C (1) 

m q 3 r 

with standard quark model parameters m q = 0.33 GeV for u and d quarks and 0.45 
GeV for s quarks, a = 0.18 GeV 2 and a s = 0.5. The wavefunctions are taken to be 
Gaussian of form exp(— p 2 /(2/3| f )) multiplied by the appropriate polynomials and 
(3m treated as the variational parameter in H for each of the IS, IP, 2S, ID states. 
Where data already exist good agreement was found. 

For transitions involving excited states with nodes in their wavefunctions there is 
a priori concern about the stability of the calculations. This does not arise in the 
transitions we consider here, and the approximation of harmonic oscillator wave 
functions is 'safe'. This can be seen, for example, in the IP — > IS transitions. 
Using the operator relation p = im q [Hr\/2, where H is the Hamiltonian (|1|), one 
obtains for the radial integrals of the wave functions the equality 

(B\p\A) = -JLu;{B\r\A), (2) 
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where A is the initial meson, B is the final meson, ui = vtia — ttib and tua and mg 
are the model, not the physical masses. With Gaussian wave functions the radial 
integrals (up to a phase) are 

(Blr\A) « (4) 

r 2 = ^ 2 + /^ 2 ) (5) 

PI = \(PI + (3 2 b) (6) 
From p[, /3 A = 0.274 GeV, (3 B = 0.313 GeV and the mass difference to = 0.562 

V ^ 7 = 0.658 (7) 

x ^ 1 7 = 0.616 (8) 

The deviation between (^) and (§) is less than 10 percent so, as the squares enter 
the widths, the uncertainty due to the use of Gaussian wave functions is about 20 
percent. This is appreciably less than the effects we are considering. 
A similar exercise for ID — » IP gives 

(B\p\A) 

V m 1 = 1.088 9 
m q 



and 



where 



and 



GeV. So we find 



and 



and 

u(B\r\A) 



1.024 (10) 



2 

so in this case the uncertainty due to the use of Gaussian wave functions is less than 
10 percent. 

The pure electric-dipole (El) transition is well-defined for heavy quarks but is cer- 
tainly a bad approximation for light quarks. We include the magnetic quadrupole 
(M2) transition as well and there is a long history of success with this approach, 
even though the M2 terms are at the same order in p 2 as El corrections proportional 
to the anomalous magnetic moment of the constituents, spin-orbit terms, Thomas 
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precession and binding effects [jH), Some of these corrections can be calculated 
14] , [X5| , p~6| , |17| |T^1 ; some can only be estimated fl4j . In any case, going beyond the 



leading approximation for the electric and magnetic amplitudes requires knowledge 
of the Lorentz nature of the confinement force. 

The most notable successes of this approach, namely of calculating the leading terms 
for each of the relevant electric and magnetic multipoles, have been in reproducing 
the magnitudes and relative phases of over 100 helicity amplitudes for photoexcita- 
tion of the proton and neutron [|I^j . These give a clear indication of which amplitudes 
are large or small, and of their relative sizes and signs. This success suggests that 
although corrections may be individually significant, their collective effect is small. 

It is possible to go beyond such an approximation, but at the price of losing some 
predictive power. Within the general assumption that electromagnetic amplitudes 
are additive in the constituents, it is possible to obtain relations among the helicity 
amplitudes, angular distributions and widths for a set of states by normalising the 
reduced amplitude to some subset of observables 
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This is the philosophy that we shall adopt. First, within the "leading multipole" 
hypothesis we can make two checks of our procedures. 

(i) We find that our result for the decay /i(1285) — > 7p is in good accord with 
experiment @, [22], ]2§. 

(ii) We predict that r(/ 2 (1270) — > 7p) < r(/i(1285) — > 7p), which appears also to 
be in accord with experiment in that there is no evidence for the radiative decay of 
the / 2 (1270) in either the MARKIII |2j or WA102 |2|, |§] experiments, and both 
have strong f 2 signals. 

Secondly, within the "single-quark-transition" hypothesis we can form a positivity 
constraint among a combination of the widths. This is satisfied by our explicit model 
dependent calculations, but enables us to draw a more general conclusion, namely 
that r(/ -> 7P) ~ r(/i -> 7p). 



Vector meson decays 

Our interest here is in the radiative decays of the p^, which we identify with the 
p(1700) 0, to / (1370) and / (1510) and of the (unobserved) (pa, to which we assign 
a mass of about 1.9 GeV, to / (1710). In the absence of glueball mixing we predict 
IH T(p D -> 7 / n "(1370)) ~ 900 keV, T(p D -> 7 /^(1500)) ~ 600 keV (assuming 



that one or the other is a pure nn state) and r(0£ — > 7/q s (1710)) ~ 200 keV. These 
widths can be changed substantially when glueball mixing is included, the degree of 
modification depending on the mass of the bare glueball. 

Three different mixing scenarios have been proposed: the bare glueball is lighter than 
the bare nn state [|25||; its mass lies between the bare nn state and the bare ss state 
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25|; or it is heavier than the bare ss state [2^]. We denote these three possibilities 



by L, M and H respectively. In |L3| we identified potentially powerful ways of 
determining glueball mixing in the scalar mesons through the radiative transitions 
Pd -> 7/o(1370) and 7/o(1500). If the bare glueball is light (~ 1300 MeV) it 
will mix strongly with the / ™ l (1370) and dilute the 900 keV width which will be 
pushed into the other scalars, in particular into the / o ' m (1500) reversing the relative 
magnitudes of the radiative widths. At the other extreme, the mixing of a heavy 
glueball (~ 1700 MeV) does not materially affect the radiative width to the / (1370) 
but severely depresses that to the / (1500) to 100 keV or less. Thus the relative 
widths of pd — > 7/0 (1370) and p£> — > 7/0 (1500) are sensitive to the glueball mass. 
There is similar sensitivity in the decays 4>d — * 7/o(1500) and 4>d — ► 7/o(1710). For 
a light glueball the former is essentially zero and the latter is predicted to be ~ 170 
keV. For a heavy glueball the situation is reversed, with the decay 4>d — * 7/o(1500) 
predicted to be about 260 keV and the decay 4>d 7/o(1710) to be essentially zero. 
These branching ratios are challenging but important as they open the possibility 
of weighing the nn and B flavour content of the scalar states and determining the 
bare mass of the scalar glueball. The only question is whether experiment will be 
sensitive to such magnitudes. 

In the isovector-scalar sector only the ao(980) is well-established. The mass, and 
even the existence, of the ao(1450) remains controversial. Its mass is a critical 
parameter in mass matrices for the mixing of the scalar glueball with the qq nonet. 
We predict T(lu d -> 7 a (1450)) ~ 610 keV and T(p D -> 7 a (1450)) ~ 85 keV. Of 
course in e + e~ annihilation or diffractive photoproduction the pn will be produced 
at approximately nine times the rate of the ljd, so the 700 (1450) rates will be the 
same for both. Observation of the ao(1450) in radiative decays would have significant 
implications, particularly when coupled with the radiative decays to the isoscalar 
scalars. 

Although our emphasis is on the radiative decays to scalars as a probe of their 
glueball content it is important to look at other decays as a check on the model. We 
find that the principal radiative decay mode of the pn is 7/i(1285), with an estimated 
width |L3[ of ~ 1100 keV. Although the decays of the /i(1285) are many-bodied, 
principally 47r and r)7nr, this is compensated for by its narrow width of 24 MeV. The 
largest radiative decay of the ujd is found to be 701 (1260), with a width ~ 1100 keV. 
Observation of these modes would confirm our calculations. The E852 experiment 
at Brookhaven has observed |27j the decay a;(1640) — > ujrj. If the ^(1640) = ujj 
then, in the 3 P model, the partial width for this mode is predicted |§ to be 13 MeV. 
Thus, if the 3 Pq model is a good guide in this measureable u(1640) — > 7a! is 

anticipated: as E852 have several thousand events in the oorj channel, it is possible 
that there could be several hundred 701 events in their experiment. 
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Scalar meson decays 



Just as the radiative decays of excited vector states to the scalars provide sensitive 
tests of their parton content, so do the radiative decays of the scalars to the ground- 
state vectors. In our model these decays are given by 



3 rriA rn 2 (3 2 ' 



r(/ - iv) = -a P —^F\i + x^-yi (11 



where A refers to the initial scalar meson, B to the final vector (V) meson, 

A = 2{i% t m (12) 
F= w exp( "5(rfbi)) (13) 

where j3 is given in fl5|) and I is an isospin factor, | for nn — > ^p and | for ss — > 70. 

From [g, p A = 0.274 GeV, (3 B = 0.313 GeV for nn states and (3 A = 0.307 
GeV, (3b = 0.355 GeV for ss states. In the absence of glueball mixing it gives 
r(/ (1370) -> 7p) ~ 2300 keV, and r(/ (1710) -> 70) ~ 870 keV. 

The width of the decay /i(1285) — > 7p is measured and provides a check on the 
model. This decay is given by 

r(A - 7 10 = |a P ^4F 2 (l + A$ + (14) 
3 mi m l p z 2 p 4 

which for an nn state predicts a width of 1400 keV. This compares well with the 
experimental value [|, E3, G3| of 1320 ± 312 keV. 



The decay ^(1270) — »■ jp is given by 

r(/ 2 _ ;V) = |ap^4^(l " 4 + (15) 
3 mi m~ p^ 10 p 4 

which is rather smaller than the others, due to the negative contribution of the p 2 / f3 2 
term. This predicts a width of 644 keV. 

Experimentally this width is small as neither the MARKIII [22] nor the WA102 
f23j experiments has any evidence for it. The branching fractions for the radiative 
decay of J/ip to /i(1285) and / 2 (1270) are comparable @, at (6.1 ±0.9) x 10" 4 and 
(1.39 ±0.14) x 10~ 3 respectively, so the non-observation of any ^(1270) signal in the 
decay J/ip — > 7(7/2) is meaningful. A similar situation holds in central production in 
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high-energy proton-proton interactions |23|] and one can deduce an upper limit 
on r(/ 2 — * jp) of 500 keV at 95% confidence level. So it is reasonable to suppose 
that our results for the fo radiative decays are valid. 

As in |H| we consider three scenarios for the glueball mass: light, medium and 
heavy. The predicted widths for the decays of /o(1370), /o(1500) and /o(1710) to 
7P and 70 are given in table 1. 





L 


p(770) 
M 


H 


L 


0(1020) 
M 


H 


/o(1370) 


443 


1121 


1540 


8 


9 


32 


/o(1500) 


2519 


1458 


476 


9 


60 


454 


/o(1710) 


42 


94 


705 


800 


718 


78 



Table 1. Effect of mixing in the scalar sector of the l 3 Po nonet for radiative decays to 
p and 0. The radiative widths, in keV, are given for three different mixing scenarios 
as described in the text: light glueball (L), medium-weight glueball (M) and heavy 
glueball (H). 

It is clear from the table that the discrimination among the different mixing scenarios 
is strong. The decay /o(1500) — > jp is perhaps the most interesting because of its 
comparatively narrow total width of ~ 120 MeV. This enhances the radiative decay 
branching fraction, which is ~ 2% for a light glueball, ~ 1% for a medium glueball 
and ~ 0.5% for a heavy glueball. The absolute magnitudes may vary by ~ 20% 
due to the uncertainties in the model, but the relative strengths manifested in the 
pattern above should be robust. Even allowing for the intrinsic uncertainties, there 
appears to be a solid conclusion that these rates are worth pursuing experimentally. 

We now assess further the robustness of these results by looking at the more general 
structures that follow from the single-quark-transition property of the dynamics. 



Single quark transition 

The results above depend upon the rather general assumption that electromagnetic 
amplitudes are additive in the constituents. Independent of details of binding dy- 
namics, it is possible to obtain relations among the helicity amplitudes, and hence 
angular distributions and widths, that depend only on the assumption that the 
mesons are described as qq states in the P and S states respectively. 

For example, in the electric-dipole approximation, the relative widths for transitions 
from fj — > 7p would be equal, apart from phase space corrections. Given the 
magnitude of one as input, the others immediately follow. While such results are a 
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good approximation for heavy flavours ||, for the lighter u, d, s flavours, significant 
magnetic-quadrupole may be expected, as is confirmed by the explicit calculations 
above. (Electric octupole contributions, while allowed in general for f 2 — > jV, will 
vanish if, as assumed here, the vector V is a pure 3 Si state, and the / 2 is a pure 3 P 2 



state PD| 



Within this assumption that the electromagnetic transition amplitude is additive 



in the constituents 0, ^TJ and that the relevant qq states can be classified as S, P 
levels, sum rules can be obtained. We can of course check that these sum rules are 
satisfied in particular explicit models. 

The most general decomposition of the helicity amplitudes for 3 Si <-> 3 Pj radiative 



transitions into electric-dipole and magnetic-quadrupole terms is given in |fL7 
The helicity amplitudes may be written as: 

3 S X <- 3 P 

A = V2{E l + 2E R ) (16) 
A = V3(E 1 + E R + M) A 1 = V3(E 1 + E R - M) (17) 



3 S 1 <-> 3 Px 



3 Si <-> 3 P 2 

A = (Pi - E R + 3M) Ai = \/3(Pi - Pr + M) A 2 = v^Pi - P fi - M) (18) 

Here Pi is the leading electric-dipole term, E R is the "extra" electric-dipole term 
and M is the magnetic-quadrupole term. In terms of our model, Pi comes from the 



convection current and E R and M from the spin-flip one (see equation (2) of [13 
The width for the decays V — > 7/j is 



r(^-7/j) = ^El^| 2 (19) 



3 A 



and for the decays fj — > is 



r(f^,V) = ^^\Aj (20) 

Consider first the decays V — >• 7/j. For equal phase space and equal form factors, we 
can eliminate the |M| 2 term by the combination 3T(V — > 7/1) — r(V — > 7/2)- This 
leaves cross terms between Pi and Pr which are in general of indeterminate sign. 
However, these too can be eliminated by including r(V — > 7/0) and a combination 
formed that is proportional to the sum |Pi| 2 + \E R \ 2 > 0. 

Thus, for equal phase space and equal form-factors, one would have 

T(p(2S) - 7/2) + 7T(p(2S) - 7/0) > 3r(p(2S) - 7/1) (21) 



8 



Similarly for the decays of pure nn fj states, 

5r(/ 2 -> 7 p) + 7T(/ - 7P) > 9r(/i -> 7 P) (22) 



It is straightforward to verify that (pip is satisfied by equations (17) to (19)0 of [T3| 
and that ( p2[) is satisfied by ([11]), (0) and (p75|) of the present paper. 



Given that empirically r(/i — > 7p) ~ 1500 keV, equation (|22f) implies that one or 
other of the f 2 or / must have a radiative width of at least ~ 1000 keV. To the 
extent that there is no clear sign of the f 2 — > 7P decay in either the data sets of 



22, 23 1 one may anticipate that the fo — > 7P could be large, in line with our specific 



calculations. In any event, such states should be sought in ip — > 77p. 



Conclusions 

Our computations show that radiative decay rates can be large for some transi- 
tions. The limitation for detection of these processes seems to be one primarily of 
acceptance. These processes should be borne in mind when detectors for future ex- 
periments, for example at Hall D in Jefferson Lab, are being designed. We advocate 
that existing data be mined to seek evidence of radiative decays, for example in 
E852 where the decay cj(1640) — > 701 should be present at a reasonable rate. If 
radiative transitions are observed, then the use of radiative decays will be proven as 
a viable technique. 

It is important to compare photoproduction and e + e~ annihilation as these give 
complementary information on the vector-meson wavef unctions. Jefferson Labora- 
tory is an obvious place for the photoproduction studies, and as regards e + e~, VEPP 
is the immediate natural candidate, though there is uncertainty as to whether its 
luminosity will be enough. DAFNE operating up to its maximum energy promises 
high luminosity in the future. Initial state radiation at BaBar, Belle and CLEO-c 
could give significant rates for e + e~ annihilation at ^/s < 2 GeV. 

A complementary approach to the role of glue in the scalar wave functions is provided 
by the decays of the scalars to 7p and 70. Both the /o(1500) and /o(1710) have 
been clearly observed in the radiative decays J/ip — > 7X, and we recommend that 
a detailed exploration of the scalar states should be planned as part of the CLEO-c 
programme. 

Radiative transitions provide a new flavour filter that can clarify the nature of the 
scalar mesons in the 1.2-1.8 GeV region, in particular the role of glue in the scalar 
wavefunctions. 



4 There is a sign error in equation (18) in [^3|. The sign of the GeGm term should be negative, 
as it is in equation (A5). 
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